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ABSTRACT 


Of the many Free Electron Lasers (FELs) in the world today, most are big and 
expensive. This is true across the wavelength spectrum, from long to very short. In 
contrast, the FEL facility in progress at NPS, which will initially operate at long 
wavelengths and at electron energies of only a few MeV, is inexpensive and smaller. 
However, longer wavelengths lead to more diffraction, which may result in beam 
spread and interaction with the undulator surfaces. Anticipating the possibility, in 
this thesis we analyze mathematically the free space Hermite-Gaussian modes of the 
optical beam, and then compare them to the Hybrid modes, where waveguide plates 
control diffraction along one axis, allowing free space diffraction along the other axis. 
We continue the analysis of the relativistic electron beam, co-propagating with the 
optical wave in the Hybrid Mode, to define new operating condition for the FEL. 
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I. 


INTRODUCTION 


The Free Electron Laser (FEL) uses a beam of relativistic electrons co-propaga- 
ting with an optical beam in a periodic magnetic field. This magnetic field causes 
the free electrons to oscillate, allowing coherent energy exchange between the electron 
beam and the optical beam. 

The first development relevant to the FEL was the microwave tube, invented in 
the 1930s, which utilized a beam of free nonrelativistic electrons in a closed microwave 
cavity. Later, atomic and molecular lasers that employed a system of bound electrons 
and an open optical resonator were devised in the 1960s. By using a beam of free 
relativistic electrons and the open optical resonator, J. M. J. Madey, the inventor of 
the FEL, combined the attributes of microwave tubes of the 1930s and the atomic 
lasers of the 1960s [1]. 

Today there are many FELs to be found in the world, and most of these pro¬ 
duce high quality laser light at IR wavelengths for scientific research. FELs generating 
short wavelengths, UV and X-ray, have been studied extensively during the past four 
decades; these require large and expensive accelerators with electron energy on the 
order of GeV. In this thesis, I will focus on long wavelength FELs, which require only 
a few MeV for the electron beam energy and much simpler technology with lower 
costs. This thesis work will provide a basic theory for designing FELs to operate at 
long wavelengths in the laboratory. Since NPS is starting a FEL facility that will 
first reach lower energies of only a few MeV, the FEL considered here could be the 
first to operate at NPS. 

The thesis is organized as follows: Chapter II describes the FEL system, briefly 
discussing how FELs work and their significant attributes. Chapter III describes the 
basic theory of FEL operation dealing with the electron behavior and interaction 
between the electron beam and the optical field. This chapter will be limited to 
plane waves to characterize the optical field. However, FELs actually operate in 
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Gaussian wave modes. So Chapter IV offers the solution of the wave equation, with 
an emphasis on Hermite Gaussian resonator modes. Then Chapter V explores FEL 
theory, including the features of Hermite Gaussian resonator modes. Chapter VI 
deals with Hybrid (HG-Waveguide) resonator modes. These are important to analyze, 
because in order to operate a FEL at long wavelengths, we may need to use waveguide 
structure capable of confining the radiation that otherwise spreads due to diffraction. 
New research is then presented in Chapter VII to suggest an FEL theory applied to 
Hybrid (HG-Waveguide) resonator modes, which is the ultimate goal of this thesis. 
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II. FREE ELECTRON LASER SYSTEM 
DESCRIPTION 


The Free Electron Laser system consists of several different elements involving 
interesting physics. Figure 1 shows the major components of an FEL system such as 
the Jefferson Lab 14kW IR FEL, and the proposed NPS FEL. There are two basic 
configurations of an FEL, the oscillator and the amplifier, which will be examined 
in the following section. Most of the components shown in Figure 1 are common to 
both designs. 

A pulsed electron beam is produced when an injector sends the beam into 
an accelerator, such as superconducting linear accelerator. The beam of electrons 
is accelerated to near the speed of light utilizing radio-frequency (RF) fields. The 
relativistic electrons enter the undulator and then get “wiggled” to produce laser light, 
bringing about coherent emission. In this process, only a small fraction (typically ~ 
1-2%) of the electron beam energy is converted to optical energy. In some cases, the 
electrons are then re-circulated through the accelerator in order to recover most of 
their remaining energy. Finally, the electron beam enters the beam dump. 



Figure 1. FEL main system schematic. Note the two configurations identified as 
the “Oscillator” or the “Amplifier.” Red dots represent the electron pulses, and blue 
ellipses describe the optical pulses. From [2]. 
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A. CONFIGURATIONS 


In an oscillator configuration, the radiation is stored and amplified between two 
reflective mirrors as shown in Figure 2. The optical beam is then outcoupled through 
a partially transmitting mirror. Coherence is developed through many passes of the 
optical field through the undulator. 



Figure 2. FEL schematic in an oscillator configuration. From [3]. 


In an amplifier configuration, the optical resonator does not exist; instead, the 
system employs a “seed laser,” which is an external laser source. Since only one pass 
is used, the amplifier uses a longer undulator to attain more gain. Coherence in this 
case is established by the seed laser. Figure 3 shows a typical amplifier configuration. 


Electron Beam 



Figure 3. FEL schematic in an amplifier configuration. From [3]. 


B. ATTRIBUTES 

The FEL has some unique attributes such as tunability, designability, high 
power, efficiency and reliability. Tunability separates the FEL from a conventional 
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laser in which the wavelength is determined by atomic energy levels, n an FEL, the 
resonance condition between the Lorentz factor 7 of the electron beam, the undu- 
lator period A 0 , and the dimensionless undulator parameter K determine the laser 
wavelength A: 

, _ A 0 (1 + K 2 ) 

X ~ ’ 

where K is proportional to the magnetic field of the undulator. This relationship 
means that by changing the spacing of the magnets, the strength of the magnetic 
field, or the energy of the electron beam, the FEL can be continuously tuned to dif¬ 
ferent wavelengths. Also, FELs have been designed over a wide range of wavelengths, 
from microwaves to X-rays; that is, they possess “designability.” The absence of a 
vulnerable medium such as gasses, fluid chemicals or solid materials means that the 
FEL has the potential to achieve extremely high power. As a result, the FEL can 
have good efficiency (~ 10%) when recirculating the electron beam. The final advan¬ 
tage of the FEL is reliability: existing FEL systems now run twenty-four hours a day 
for weeks. On the down side, the FEL is as huge as a laboratory and considerably 
expensive. 

C. COMPONENTS 

The major components of an FEL system are an injector, an accelerator, an 
undulator, a resonator for the oscillator, a seed laser for the amplifier, and a beam 
dump. There are many other constituents that allow the FEL to operate, such as 
bending magnets for steering the path of the electron beam, an optical beam transport 
system, radiation shielding and vibration control. 

1. Injector 

The injector generates and accelerates the beam of free electrons to about ~ 5 
MeV energy. There are several ways to produce the beam of free electrons. Typically, 
the free electrons are emitted into vacuum through either thermionic emission or the 
photoelectric effect. In the thermionic cathode, the electrons are emitted from the 
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surface of the cathode when the cathode is heated sufficiently. A photocathode uses 
an incident drive laser to excite the electrons in a cathode. By using the photoelectric 
effect, electrons are ejected as long as the photon energy is greater than the cathode 
work function. The released electrons can be emitted into a Radio Frequency (RF) 
cavity, that accelerates the electrons to the desired energies (~ 5 MeV) before they 
enter the main accelerator. The electron beams usually come out in pulses with 
excellent quality (a narrow, collimated beam with low energy spread); not maintaining 
such beam quality can affect FEL system performance. An RF injector is shown in 
Figure 4. 



Figure 4. Schematic of an injector. From [4], 


2. Accelerator 


After the injector, the electron beam is now passed on to the accelerator. 
The accelerator increases the energy of the electron beam from ~ 5 MeV to the 
desired level. Direct Current (DC) acceleration is possible, but most FELs use a 
Radio Frequency Linear Accelerator (RF LINAC) to obtain a final energy of ~ 100 
MeV. The accelerator may consist of several evacuated superconducting metal RF 
cavities. The electromagnetic field generated in each cavity should be in phase with 
the injected electron pulses to provide net positive acceleration. 
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3. 


Undulator 


The undulator is a significant part of the FEL located inside the optical res¬ 
onator cavity and typically consists of alternating permanent magnets. When the 
electrons leave the accelerator, they enter the periodic magnetic field and are made 
to wiggle in a transverse direction by the Lorentz force. The wiggling motion causes 
the electrons to radiate photons and amplify the fight already in the undulator. As 
the electrons and light interact in the undulator, there is an exchange of energy: some 
electrons gain energy from the fight and some lose energy to the fight. The electrons 
that gain energy move a little bit faster along the axis of the undulator and those that 
lose energy move a little bit slower; this causes “micro bunching” of the electrons on 
the scale of an optical wavelength. Undulators typically have either linear or helical 
polarization. A linear undulator drives the electrons to move sinusoidally in a plane 
while, a helical undulator drives them to a corkscrew motion. A schematic of an 
undulator and the electron beam path is shown in Figure 5. 



Figure 5. Schematic of an undulator. The red line is the electron beam wiggling along 
the transverse axis. The yellow fines indicate the stimulated emission propagating 
along an axis of forward direction. The undulator consists of magnets alternating 
north and south poles, resulting in the periodic magnetic fields. From [5]. 
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4. 


Resonator 


In the oscillator configuration, the resonator is the component in which the 
light is stored and amplified over many passes. It is made up of an evacuated cavity 
bounded by two mirrors. One of the mirrors has high reflectivity, while the other 
allows a percentage to come out through a partially-transmitting mirror. The light 
derived from this out-coupling mirror is the usable laser light of the FEL system. It is 
obvious that for maximum coupling efficiency, electron pulses must be synchronized 
with the bounce time of optical pulses between resonator mirrors. A schematic of the 
electron and optical pulses passing through the undulator in the resonator is shown 
in Figure 6. 


electron beam pulses 

wiggler 

f -mm* 


100 % 


re£1acting 
mirror 


• - 

* 

partially 

transmitting 

mirror 


optical 

laser 

pulses 


Figure 6. Schematic of a resonator. From [2]. 


5. Beam Dump 

The beam dump is the last part of the electron beam path. It consists of an 
absorbing material (such as a large metal block) surrounded by radiation shielding. 
Extracting energy from the electron beam amplifies the optical beam, but the extrac¬ 
tion of energy is typically only a few percent, leaving a high-energy electron beam 
which results in poor system efficiency, and lots of potentially dangerous radiation 
emitted at the dump. If the system uses recirculation of the electron beam, the elec¬ 
trons enter to the LINAC a second time 180° out of phase, decelerating electrons, 
and giving most of their energy back to the RF field. The recirculated electron beam 
has approximately ^3 MeV average energy at the end of the LINAC, having lost a 


few percent to the laser beam and leaving about the same energy that the electron 
beam had when it was generated at the injector. This greatly improves overall system 
efficiency, and reduces radiation shielding requirements at the dump. A schematic of 
a beam dump is shown in Figure 7. 



Fe (2-3m) 


Cu (5m) 


Surface 

91°C 


Seam 


Taper cut 360°C 
((|)40cm x 4m) 


damp250-52 Damp(Cu-Fe) @ 3.G&l4ppp & 6D0W/m£/K tp ^Ocm-Gaussiar. 


Figure 7. Schematic of a beam dump. From [7]. 
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III. 


BASIC FREE ELECTRON LASER THEORY 


We will now explore the mathematics that describe the operation of the FEL. 
The physics of the electron’s behavior and the interaction between the electron beam 
and optical field are pivotal to the operation of FEL systems. In the following sec¬ 
tions, the x and ?/-axis designate the transverse axes while the axis designates the 
longitudinal axis passing through the center of the undulator. 

A. ELECTRON BEHAVIOR 

The electrons in the undulator have a microscopic evolution. 

1. Undulator Field and Transverse Motion 

Exploring the behavior of electrons in the magnetic field of the undulator is 
necessary to understand the physics of the FEL interaction. We will deal with rela¬ 
tivistic electrons entering a helical undulator field (which is simpler mathematically 
than a linear undulator) in helical orbits around the 2 -axis in the presence of the 
optical field. An electron passing through the undulator is exposed to three fields: 
the magnetic field from the undulator magnets, the magnetic field of the optical wave 
and the electric field of the optical wave, as shown in Figure 8. 



Figure 8. For convenience, undulator and Optical Fields in the case of linear undula¬ 
tor. From [2] 
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The magnetic field of the helical undulator is 


B u = B [cos (koz ), sin (koz ), 0] 


tun) 


where B is the magnetic field amplitude, k 0 = 2 tt /A 0 is the undulator wave number 
and A 0 is the undulator’s period which is the distance of one complete cycle of the 
magnetic field in the undulator. The corresponding electric and magnetic field of the 
optical field are 

E r = E (cos ip, — sin ip, 0) (III.2) 

B r = E (sin ip, cos ip, 0) (III.3) 


where E is the optical wave amplitude in cgs units, ip = kz — wt + (p, k = 2ir/X is the 
optical wave number, w = kc is the optical angular frequency, and <p is the optical 
phase. 

The relativistic Lorentz force equations determine the electron motion in the 
undulator. The relativistic Lorentz force equations in cgs units [1] for an electron are 


^ (7/3) = [E r + /3 x (B u + B r )], 
dr me 


d7 
d t 


e 

me 


P • E r , 


7 2 = 1 -/ 3 -/ 3 , 


(111.4) 

(111.5) 

(111.6) 


where v = /3c is the electron velocity, e =| e | is the electron change magnitude, m is 
the electron mass, c is the speed of the light, and 7 is the relativistic Lorentz factor. 
Substituting the fields from Equations (III. 1 ), (III. 2 ), and (III.3) into the (III.4), we 
get 

-77 (7/3®) = —— [E cos ip (1 - /3 Z ) - (3 Z B sin (k 0 z)\ , (HI-7) 

dr me 

-77 (7A/) = —““ [-Esin^ (1 - f3 z ) + (3 Z B cos (k 0 z)] , (III.8) 

ar me 

d 6 

—r, (7 Pz) = - [E (f3 x cos ip - (3 y sin ip) + B{f3 x sin (k 0 z) - f3 y cos (k 0 z)}]. (III.9) 

ar me 
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Since the electrons are traveling close to the speed of light along the z axis, /3 Z ~ 1, 
so that terms containing (1 — (3 z ) can be ignored in the transverse x-y components of 
Equations (III.7) and (III. 8 ). In other words, E (1 — /3 Z ) is small compared to f3 z B, 
which means magnetic fields derived from an undulator have a bigger effect on the 
electron’s transverse motion than the optical fields. Integrating Equations (HI.7) and 
(III. 8 ) with respect to time and setting the constant of integration to zero, indicating 
perfect injection into helical orbits, we get 


— cos (koz ), 

7 

(m.io) 

— sin ( k 0 z ), 

7 

(uni) 


where K = eB rrns \ 0 /2irmc 2 is the dimensionless undulator parameter and B rms is 
the rms value of the undulator field. For the helical undulator B rms = B] for the 
linearly polarized undulator B rms = B/\[2. Integrating once again and using k^z = 
kov z t = ko(3 z ct ~ k 0 ct = cco t, we get 

rK K An 

x ( t ) =-sin ( u> 0 t ) =-sin (cco t ), (III. 12) 

cuo7 27T7 

cK K An 

y ( t ) =-cos (uot) = -cos (coot ), (III. 13) 

w'n'; 27T7 

which describes the transverse (helical) motion of the electron. 


2. Longitudinal Motion and the Electron Phase 

Now consider the effects of the optical electric field. We scrutinize the micro¬ 
scopic electron motion along the longitudinal axis which, while small in magnitude, 
is key to FEL operation. Substituting Equations (III.10) and (III. 11 ) into (III.9) and 
using a trigonometric identity gives 

7 Pz = —— cos ( k 0 z 4- ip) - 7 / 3 Z . (III. 14) 

7771 c 

In order to get 7 , substitute Equations (III.10), (III. 11 ) and (III. 2 ) into (III.5) and 
get 


dq 

df 


= 7 


eKE 
7 me 


cos (koz + 'ip). 
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(III.15) 



Substituting Equation (III. 15) into (III. 14), we arrive at 


• ( r\ r, 

f3 z = (l-f3 z )— — cos (k 0 z + *ip). (III.16) 

7 z mc 

Also, we can rewrite the Equation (III. 15) by using the definition of pj and the “elec¬ 
tron phase,” C = (k + k 0 ) z — out, which describes the electron’s phase with respect to 
the combined optical and undulator fields, to obtain 

_ eKE_ cog (^ _)_ 0) . (III.17) 

qmc 

Since ho, k and ou are fixed, the dynamical variable in ( (t) is 2 (t) which represents 
the electron’s microscopic position, ( cx kz (t). Note that when cos (( + (/)) > 0, the 
electron energy increases by absorption, 7 > 0; when cos ((' + 7 ) < 0, the electron 
energy decreases by stimulated emission, 7 < 0. This causes the electron beam to 
bunch and to radiate coherently. 

3. Resonance Condition and Pendulum Equation 

FEL resonance occurs when the energy exchange between the electron and the 
optical beam is optimum. The electron travels down the z-axis with a speed (3 z c while 
the optical pulse travels at the speed of light, c. The relationship between these two 
speeds brings about the “electron-photon race,” which is illustrated in Figure 9. 



^0 ^ 


Figure 9. Electron-photon race. From [2]. 

The time for an electron to travel through one undulator period is Ao / /3 z c. During that 
time, one wavelength of light is emitted because the electron executes one oscillation 
in the undulator. The wavelength of light then travels a distance A + A 0 at speed c. 
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So, 


Ap 

P Z C 


A + Aq 


From (III. 18) we see that the laser wavelength at resonance is 

Ao (1 — Pz) 


A = 


Pz 


(III.18) 


(III.19) 


Now we hnd an expression for (5 Z in terms of the dimensionless undulator 
parameter, K. From Equations (III. 10) and (III. 11 ), we find that /3j _ 2 = K 2 /q 2 . 
Therefore, the relativistic Lorentz factor, 7 , which indicates the relationship between 
the electron energy qmc 2 and the z-axis velocity c/3 z can be expressed as 


K 2 


7 “" = 1 - P = 1 - P± ~ P z = 1 - — - Pz- 

7 2 

For relativistic electrons, 7 1, we can use the binomial expansion to get 

1 + K 2 


(III.20) 


Pz 


1 - 


2 7 2 


(III.21) 


We can rewrite the resonance condition, Equation (III.9), in terms of K to get 

A 0 (1 + K 2 ) 


A 


2q 2 


(III.22) 


By changing the physical parameters, such as A 0 , K and 7 , we can get a wide range of 
wavelengths of the optical beam, which indicates the FEL’s tunability. Also we can 
now see that the wavelength of the optical beam A is much smaller than the undulator 
period A 0 because of the small factor I/ 7 2 . 

Now we are ready to derive the FEL “pendulum equation” by combining some 

o 

of the previous equations. Let us introduce an operator (..)= d(..)/ dr which can 
be expressed as 


_d_ 

dr 


d d t 


L _d 
c dt’ 


dt dr 

where t = ct/L is the “dimensionless time,” where L is the undulator length, so that 
r = 0 —> 1 along the undulator. By using this relation, we can determine the rate of 
change of the electron phase along the undulator as the electron phase velocity, 

dC 


v — ——C— L [(k + ko) P z — k ], 


(III.23) 
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where the derivative of ((electron phase) is with respect to r; again L = NX 0 is the 
length of the undulator, so that 0 < r < 1. By substituting Equation (III.21) into 
(III.23) and using k k 0 near resonance, we arrive at [6] 

dv = 4yriV—, (III.24) 

7 

relating changes in beam energy to the electron phase velocity. Combining Equation 

(III. 15) and (III.24) gives 

dv L dv L dv dq LAirNeKE 

—— =-— =---— =-cos 

dr c dt c dq df c q qmc 

v = C= | a, | cos (C + <j>) j (III.25) 

where \ a \ — dirNeKLE/^mc 2 is the “dimensionless optical field amplitude.” Equa¬ 
tion (III.25) has the form of the pendulum equation, and describes the microscopic 
motion of the electrons in phase space, (£, v). 

4. Phase Space Evolution 

Since phase velocity is associated with the electron energy according to Equa¬ 
tion (III.24), phase space diagrams show an energy exchange between electrons and 
the optical field. Assuming | a \ and q are approximately constant, the electrons 
follow phase space paths given by [8] 

v 2 = vq + 2 | a | [sin (C — o) — sin (( 0 + 4>)}, (III.26) 

where (, 0 and z/ 0 are the initial electron coordinates in phase space. A special phase 
space path is the FEL separatrix [9], starting at the unstable fixed point (£ 0 = dn/2 
and vo = 0), so that 

v 2 = 2 | a | [1 + sin (£ + 0)]. (III.27) 

From this equation, we can see that the peak to peak height of the separatrix is 
41 a I 1 / 2 . A larger | a | results in the phase space paths evolving much faster and a 
separatrix increased in height. Inside the separatrix, the electrons follow closed orbits 
(corresponding to a pendulum swinging back and forth). Outside the separatrix, the 


(C + 0 ) > 
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electrons follow open orbits (like a pendulum swinging over the top). Figure 10 shows 
the evolution of 20 sample electrons at resonance z/ 0 = 0 and weak field a 0 < 7r during 
the dimensionless time from r = 0 (yellow dots) to r = l(red dots). 

0.2 


0 

0.1 

0 

-1t/2 £ 3JC/2 0 X 1 

Figure 10. Phase space evolution of electrons injected at resonance (uq = 0), and 
corresponding optical field gain and phase evolution. From [4]. 


0.2 

0 

0.1 

0 

-ji/2 C , 3rc/2 0 t 1 

Figure 11. Phase space evolution of electrons injected off resonance (uq = 2.6), and 
corresponding optical field gain and phase evolution. Periodic boundary conditions 
are applied to keep the electron phase £ between — 7 t/2 and 37r/2. From [4], 




The symmetrical distributions of the electrons indicates that about half of the elec¬ 
trons gain energy from the optical field (paths shift up) and, the other half lose energy 


17 





























to the optical field (paths shift down). We can see a bunching of electrons occurs but 
there is no gain since an equal number of electrons gained and lost energy. 

If electrons were injected slightly off resonance (uq = 2.6), as shown in Figure 
11, more electrons would shift to a lower value of phase velocity than shift to higher 
energy. There is an obvious bunching of electrons near ( = n, which causes the 
positive gain. 

Gain is the ratio of the change to the initial value in the power of the optical 
field as it passes along the undulator. We can express the gain as 

a 2 (r) — ao 2 


G t) = 


a 0 Z 


(III.28) 


where a(r) is the field along the undulator at r and ao is the initial field starting at 
beginning of the undulator. 


B. ELECTRON-OPTICAL FIELD INTERACTION 

We derived the pendulum equation which represents the microscopic evolution 
of the electrons in the FEL interaction region. We now need to derive an equation 
describing the self-consistently evolving optical field. 

1. Maxwell’s Equations 

Now, consider the interaction between the electron beam and optical field, and 
start with Maxwell’s equations in cgs units [10]: 

V • E = 4-7rp 


VB = 0 

„ „ 13B 

V x E = — — 
c dt 


4tt t 1 <9B 
V x B = —J + 

c c at 


(III.29) 

where p is the charge density, J is the current density. From Maxwell’s equations, we 
can derive the full wave equation associated with transverse current [8], 

1 cP“ \ Air 

* 2 ~*w) am = -T 3±m ' (IIL30) 
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where A(x, t) is the optical vector potential, c is the speed of light, and J_l(x, t) is the 
transverse current density; in our case, due to the wiggling motion of the electrons 
traveling through the undulator. 

2. FEL Wave Equation 

For our helical undulator, the optical vector potential can be written as 

A(x, t) = ie ia , (III.31) 

k 


where E =| E | is the complex laser electric field, a = kz — cut is the carrier 

wave, £ = (—*,1,0) is the laser field’s polarization vector for the case of circular 
polarization from a helical undulator, and x = xi + yj + zk is the position vector. 
By substituting Equation (III. 31) into (III. 30) and assuming that the optical wave 
amplitude | E(x, t) | and phase </>(x, t) are slowly varying along the z axis, we arrive 
at [2] 


o i d Id 

V i + 2 ik ( — I — 
oz c at 


E = J ± -£*e~ ia , 


(III.32) 


which is the parabolic or paraxial wave equation with a source current Jj_. When we 
change the variables to a reference frame that “follows the light,” u = z — Lt = z — ct 
and r = ct/L, then we can rewrite the operator [11] 


d Id Id 
dz cdt L dr 


(III.33) 


The source current which is the sum of all single-particle currents can be expressed 
as 

J_i_ = —ec ± h (3) (x - x;(f)), (III.34) 

i 

where h®(x — x*(t)) = S(x)S(y)S(z) is the three-dimensional Dirac delta function, 
and Xj(£) is the i th electron position at time t. The transverse motion of an electron 
(3 ± which contributes to the transverse current in the undulator can be expressed 
as 

f3 1_ = —— (cos k 0 z, sin k 0 z, 0) = Re ( —— i£e~ lkoZ ) . (III.35) 

7 \ 7 / 
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Substituting Equation (III. 35) into (III. 34), and (III. 34) into (III. 32) gives 


V ± 2 + 


2 ik d 
~Tdr 


E(r,r ) = —47 rieKkp(x.,t) (e ^/q) , 


(III.36) 


where r = xi+yj , and p(x, t) = f 8^ (x — Xj (t)) dV is the local electron density in 
a small volume element dV, and < ... > indicates an average over sample electrons in 
the volume element. If we multiply Equation (III. 36) by (AirNeKL/^ 2 mc 2 )(—Li/2k) 
and define the complex laser field a =\ a \ with the previous laser field amplitude 


a \ = Air N eK LE / ^/ 2 mc 2 , we can get 


iL„ 2 d 

2k 1 dr 


a 


(r,r) = -(je * c > , 


(III.37) 


where j = 8 N (e tt KL ) 2 p/ y 3 m c 2 is the “dimensionless current density” which repre¬ 
sents another important FEL parameter (maybe the most important). By introduc¬ 
ing the dimensionless transverse coordinates x = x (k/2L ) 1 ' 2 , y = y(k/2L) 1 ' 2 and 
replacing the Laplacian operator with a dimensionless Laplacian operator, Vj_ 2 = 
(k/2L)\7 2 ± , we arrive at the FEL wave equation [11] 

d 


:Vi + 


dr 


t (r,r) = - (je * c ) , 


(III.38) 


where = d + dy 2 , and the coordinates ( x,y,z,r ) are all now dimensionless. 
When diffraction is neglected, the FEL wave equation can be written in the simplest 
form 

°a= -j (e“ iC ) . (III.39) 


The dimensionless current density j measures the coupling between the electron beam 
and the optical mode, and (e - ^) measures the amount of electron bunching. When j 
is large (j 3> 7r), the coupling is large and we have high gain; when j is small).) < 7r) 
we have low gain. 


3. Gain 

Until now, we have studied the interaction between free electrons in the beam 
and laser light through the Lorentz force equations in the helical undulator. Previ- 
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ously, we have understood the bunching of electrons in phase space resulting in coher¬ 
ent radiation through the pendulum equation. Now, we have coupled this bunching 
to the laser beam through the FEL wave equation. 

When we assume a weak field | a, |< 7r and low gain j < tt, we can apply a 
pertubation expansion to the pendulum equation, resulting in FEL gain, [2] 

2 — 2cOs(zqr) — U 0 TSm(u 0 T ) 

^o 3 

where a filling factor F represents a ratio of cross sectional areas of electron beam to 
optical mode such as F = (ry/wq) 2 at the waist of the mode, where ry, is the electron 
beam radius and w$ is the optical mode waist radius. Gain is proportional to the 
dimensionless current density j and filling factor F, and the final gain spectrum at 
r = 1 depends only on the initial electron phase velocity zq as shown in Figure 12. 
As seen before, there is no gain at resonance, zq = 0, and the peak of G/jF (gain 
normalized to j and F) is about 13.5% and occurs at zq ~ 2.6. // 


(III.40) 


G{r) = 


a ( T ) 

a 0 2 


- 1 = jF 



Figure 12. Single pass G/jF versus initial electron phase velocity zq at r = 1. 
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4. Extraction 

In order to describe the energy exchange between the electron beam and optical 
field, we estimate the extraction r\ defined as the fractional energy transferred from 
the electron beam to the optical field during a single pass through the undulator [6]: 


V = 


_ Eppt 


(A E) (Ay) 


(III.41) 


E b E b 7 

where E opt and E b are the energy in the optical field and the electron beam, respec¬ 
tively, in a common volume element of the beams, and A E = Ayrric 2 is the net 
change in the electron beam energy through the undulator. Applying Equation (III. 
24) to (III. 41) allows us to have the following result: 


_ (Ay) _ (A u) _ J_ 

Vmax y 4 t tN ~ 2 N 


(III.42) 


where it has been assumed that approximately half the electrons lose energy across 
the positive gain bandwidth Av ~ 7r illustrated in Figure 12. We know the height of 
the separatrix is 4 y/a, as presented below Equation (III. 27). This height corresponds 
to a change of an electron phase velocity Au, when an electron is in the largest closed 
orbit. Suppose that only half of the electrons are in closed orbits, then (Av) would be 
2\fci. Now, we need to prove that a ~ tt 2 by giving a short account of “Synchrotron” 
oscillation. 

Another argument starts from phase space considerations. Start with the 
pendulum equation with low gain j < n and strong fields | a |> tt. Assuming | a |~ 
constant and 0 = 0 yields the pendulum equation as 


C= I a I cos (C) • 


(III.43) 


For an electron trapped in strong field near the stable fixed point (Co = "/2, 0) in the 
phase space, its motion can be described by 


CO) = ! + 


(III.44) 
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where Co = tt/2 is the initial electron position and x(t) <C ir [2], Substituting Equa¬ 
tion (III. 44) into (III. 43) gives 


OO | 

x=\ a 



a x = —v o x , 


(III.45) 


where v s is the synchrotron frequency about the stable fixed point £ 0 = 7t/2. We 
recognize that the simple harmonic oscillation of electrons around the stable fixed 
point (7t/2,0) has the frequency u s = \ a | 1//2 . The trapped electrons will execute 
approximately half a synchrotron oscillation at saturation, or v s ~ 7r when | a |~ ir 2 . 
Assuming about half the electrons are trapped, we arrive at the estimate for extraction 


= ( A 7) = (A^) _ Va _ J_ 
7 4^^ 47rAi 2 N 


(III.46) 
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IV. 


HERMITE GAUSSIAN RESONATOR MODES 


In the previous chapter, we examined the interaction between the electron 
beam and the optical field inside the undulator. This chapter examines the wave 
nature of laser beams as they propagate in free space. Up to this point, we have 
assumed the optical beam to consist of plane waves neglecting diffraction. The wave 
usually takes the form of a Gaussian beam as defined by the optical cavity bounded by 
spherical resonator mirrors. In order to understand the characteristics of a Gaussian 
beam, we will find solutions to the free-space wave equation including diffraction. 


A. SOLUTION OF WAVE EQUATION 

The characteristics of Gaussian laser beams in a resonator have been analyzed 
by Siegman [12], Kogelnik [13], Pedrotti [14], and others. This chapter follows their 
analysis. 


1. Scalar Wave Equation 

Let us consider a general wave equation for the vector potential A without a 
current source, 


V 2 -^)a = 0 or VA=I^. (IV.l) 

Assuming A(x,y,z,t) = u(x,y, z)T(t) and utilizing a process of separation of vari¬ 
ables, we have 


1 11 d 2 T 

—V 2 w = — — —— = constant = — k 2 , k > 0, (IV.2) 

ii. c 2 T dbg 

space time 

where the left-hand side is only a function of space and the right side is only a function 
of time, so that u represents a component of the electromagnetic field or of the vector 
potential A. When we consider only the scalar function of space, u, we arrive at a 
scalar wave equation, the Helmholtz equation [15] 


(V 2 + fc 2 ) u(x, y, z) = 0, 


(IV.3) 
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where k = 27t/A is the propagation constant in the medium, with a harmonic time 
variation, e~ luJt . 

2. Fundamental Mode 

For an optical beam traveling along the ^-direction, Equation (IV.3) becomes 

u(x, y, z) = p(x, y, z)e lkz , (IV.4) 


where p is a slowly varying complex valued function, and e lkz describes a plane- 
wave in the ^-direction. Generally, p represents the distinction between the actual 
laser beam and a plane wave, which we will examine in detail below. Substituting 
Equation (IV.4) into the scalar wave Equation (IV.3) yields, in cartesian coordinates, 
the reduced equation 


d 2 p d 2 p d 2 p> , dp 

d^ + w + w + & =a 


(IV.5) 


As z ranges over an optical wavelength the change in p will be small, also the variation 


of p with respect to z will be much smaller than with respect to x and y. This 
slowly varying dependence of p on z can be expressed mathematically by the paraxial 
approximation [12]: 


d 2 xjj 

<c 

2 


d 2 ^ 

, and 

d 2 ^ 

dz 2 


dz 

5 

dx 2 


dy 2 


(IV.6) 

We can then discard the second derivative of p with respect to z in Equation (IV.5) 
and obtain the paraxial wave equation: 


d 2 p d 2 p dp 
dx 2 dy 2 + dz 


(IV.7) 


A form of solution for Equation (IV.7) in the fundamental mode can be written as 
[13] 


p(x,y,z) = exp 


P(z) + 


k 


2 q(z) 


(IV.8) 


where r 2 = x 2 + y 2 , P(z) is a complex phase shift which is related to the propagation 
of the optical beam, and q(z) is the complex curvature of the phase front which 
represents the spherical shape of the Gaussian beam. We need to find these unknown 
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functions, P(z) and q(z). In order to get each term of Equation (IV.7), we take 
derivatives of p) with respect to each component. We then have the following sources 


d 2 p) 

dx 2 

d 2 p) 

w 


ik k 2 9 . , 

--z 2 P), 

q q z 

ik k 2 9 . , 

- -y P), 

q q z 


(IV.9) 
(IV. 10) 


and 


dp) _ f.dP 

dz \ dz 


ik 


p). (IV.ll) 

Inserting Equations (IV.9), (IV.10), and (IV.ll) into (IV.7) and dividing through by 


_ 

2 q 2T dz 


p), we arrive at 


A r 2 (<% _ ! 

q 2 \dz 


\dz q 


(IV. 12) 


In order for Equation (IV.12) to be valid for any r 2 , the expressions inside both 
parentheses must be zero. Therefore, we obtain the relations 


dq 

dz 


(IV. 13) 


and 


(IV. 14) 


dP _ i 
dz q 

At this point, we should introduce two important parameters of the Gaussian 
beam, w(z) and R(z), where w(z) is the beam radius or spot size, and R(z) is the 
radius of curvature of the wavefront as shown in Figure 13. These parameters are 
associated with the waist spot size, w 0 , and the ratio zjz 0 by the formulas [12] 

.2 


w(zY 


w 0 


1 + ^2 
Zo 2 


(IV.15) 


which can be derived from Equation (IV. 13) and 


R{z) — z ( 1 + 


A) 


(IV. 16) 


where zq = ttu ! o 2 /A, the Rayleigh length, is the distance which the beam travels from 
the waist to where the beam area doubles. 
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Figure 13. Notation for a fundamental Gaussian beam diverging away from its waist. 
w 0 is the waist spot size, w(z) is the beam radius, and R(z) is the radius of curvature. 
After [12] 

We now consider Equation (IV. 13); integrating this equation gives 

q = z + q 0 = z - iz 0 , (IV.17) 

where the initial value is qo = —mwo 2 /\ [13]. Starting with Equation (IV.17) and 
using Equations (IV. 15) and (IV. 16), we get the following result: 


1 

q 


z 2 + z 0 2 


+ i- 


Zq 


+ z 0 2 


1 A 
= -^ + i -; 


R 


7 TW* 


(IV. 18) 


(IV. 19) 


For convenience of later calculation, we can rewrite this equation as 

k k ikX k i 

2 q 2 R 2itw 2 2 R w 2 ’ 

where k = 2tt/A. To determine P(z), we insert Equation (IV.17) into (IV.14) and we 

have 

dP iz z 0 

dz Z 2 + Zq 2 Z 2 + Zq 2 ' 


(IV.20) 
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Integration of this equation yields the following result 


P{z ) = i 



-i ^ 

tan —. 

z 0 


By using Equation (IV. 15) we can rewrite this as 


(IV.21) 


iP(z) = In (—) - itan- 1 —, (IV.22) 

V w / z o 

where the tan -1 (z/zq) represents the phase shift difference between the Gaussian 
beam and a plane wave, and the amplitude factor, vjq/vj, describes the decrease of 
the beam amplitude on the axis because of the expansion of the beam. Substituting 
Equations (IV. 19) and (IV.22) into the trial solution, Equation (IV.8), we can finally 
obtain the solution of the fundamental mode: 


w 


exp 


kr 2 
~2R 


-i 


z o 


(IV.23) 


Using this solution, Eqnation(IV.4) can be written in the form 

, < ' 

, x Wq 2/ 2 

u(r,z) = —e ' exp 


w 

Amplitude 


kz 


kr 2 


— tan 1 — 
zo 


z 


2 R 


V 


-V- 

Phase 


1 


(IV.24) 


Gaussian beams do not change their profile, but beams which have long wave¬ 
lengths will be expanded significantly due to diffraction as they propagate. Figure 14 
shows well the expansion of Gaussian beams. In order to show the expansion, we set 
the wavelength to A = 300pm. 


3. Higher Order Modes 

In the previous section, we discussed only the lowest-order solution of Equa¬ 
tion (IV.7). There are other higher-order solutions for the Equation (IV.7); these 
can take the form either of Hermite-Gaussian functions in cartesian coordinates or 
of Laguerre-Gaussian functions in cylindrical coordinates [12]. The properties of 
Laguerre-Gaussian modes in the FEL were analyzed by A. Kampouridis [4]. Through¬ 
out the rest of this thesis, we will use the cartesian coordinates for consistency. 
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Figure 14. Surface plots of a Gaussian beam evolution as propagating along the z- 
axis. When we assume 2=lm, (a) is an intensity profile at the beam waist (z=0.5m). 
At (b) (z= 0.75m) and (c) (z= lm), we can see the expansion of the beam as it passes 
along the 2-axis from the waist. 
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Let us assume a trial solution for Equation (IV.7) as 


/ ( x \u(V\ T- {nr \ , k 2 
V = Q — )h - exp i\ P(z) + ——r 

\ws \wJ l \ 2 q[z) 


(IV.25) 


where <7 is a function of x and 2 (because w depends on z), and h is a function of 
y and z to account for diffraction. Now, we should find g(x/w), h(y/w), P(z), and 
q(z) to get the complete solutions. 

For convenience, let us find derivatives of g and h with respect to x and y. By 
using the chain rule, we get 


dg dg d{x/w) 1 dg 1 , 
dx d(x/w) dx wd(x/w) ’ 

d 2 g d fdg\ 1 fdg'\ 1 dg' d(x/w) 1 


dx 2 dx \dx J w\dx) wd(x/w) dx w 2 ^ 


(IV.26) 


1 \ dw x ,dw 


dg dg d(x/w ) dg / 1 \ dw x ,dw 

dz d(x/w ) dz X d(x/w ) \ w 2 ) dz w 2 ^ dz ’ 


--7i . 7W = ,>'/ ■ (IV.27) 


where (..)' = d(..)/d (x/w). We can easily find the derivatives of h through similar 


processes: 


dh 1 , d 2 h 1 „ 

dy w ’ dy 2 w 2 

dh y , dw 

dz w 2 dz 


(IV.28) 


(IV.29) 


We are ready to take derivatives of Equation (IV.25). After some tedious work, we 


obtain the following results: 


d PL = e (<r + ^) J_ g% + l g , h (+ gh hi) _ gh (t x 


w \q 


(IV.30) 


V' + V ( -y) +#(-)- gh( k -y 

dy z w z w \ q J \ Q J \Q 


(IV.31) 


£ = eW) - lyh'^y + . (IV.32) 

dz \ w 2 dz w 2 dz dz 2 q z dz J 


31 



Substituting Equations (IV.30), (IV.31), and (IV.32) into (IV. 7) and dividing by 
gh exp (iP + T^r 2 j, we arrive at the following complicated equation 


1 g" 2 ikx g' 


w z g' 1 h" 2 iky h! 


w z h! 
w 2 h 


w 2 g wq g w 2 g w 2 h wq h w 2 

'-v-' '-v- 

independent of y = — u x independent of x = —u 2 


r 2 ( — 


k 2 k 2 dq 
q 2 q 2 dz 


r 2 only = 0 


2 ik dP 

- 2k ir 

q oz 

'-V-' 

no r = «i + m 2 


(IV.33) 


In order for this equation to have a solution, let us assume that the ^-independent 
terms are —U\ and the ^-independent terms are — u 2 . The r 2 -dependent term should 
be identically zero since it must be valid for an any location in the transverse plane 
[11]. In the end, the terms left at the end should be u\ + u 2 so that the left-hand side 
sums to zero. 

Let us consider the r 2 terms of Equation (IV.33) to get 


, k 2 f dq 


q 2 \dz 


1 = 0 . 


(IV.34) 


We then have the familiar relation, ( dq/dz ) = 1, which we have seen in the funda¬ 


mental mode derivation. Thus, we have, easily: 

1 _ 1 .A 
q R+ irw 2 


(IV.35) 


Next, multiplying the ^-dependent terms of Equation (IV.33) by gw 2 and 
collecting terms in the derivatives of g, we have 


// ^ , ( w dw \ . o 

g + 2 ikx - — ) g + u x w 2 g = 0. 

\q oz J 


(IV.36) 


Let us calculate only the 2nd term. Recalling Equations (IV. 15) and (IV. 16) and 
substituting these two relations into the 2nd term, we have 


, (w iX dw\ , 2-7T / iX 

2 ikx — H- — g = 2 ix— — 

V K 7 TW OZ / A V 7TW 


32 



where dw/dz = w/R. We can easily figure out this relation when we take the 
derivative of w with respect to We can then rewrite Equation (IV.36) as 


9 


// 


^±.1 . 2 „ 
—g + utw g = 0, 
w 


(IV.37) 


which is close to a form of the Hermite equation y" — 2 xy' + 2 my = 0 [16]. In order 
to transform Equation (IV.37) into the precise form of a Hermite equation, let 


X = 




(IV.38) 


Then we have 



dg 

d (x/w) 



Inserting Equation (IV.39) into (IV.37) yields 



d 2 g 

d X 2 



+ 2mg 


0, 


(IV.39) 


(IV.40) 


which is exactly the same as the Hermite equation where 2m = u\w 2 /2. If m = 
0,1, 2,..., then a solution of the Hermite equation is the Hermite polynomial H m (x) 
given by [16] 


six) = H m (x) = H m (^) = (-ire 2 ^; (e-^) • (IV.41) 

There are similar solutions for the y dependent terms of Equation (IV.33), i.e., h(£) = 
H n (£) = H n (y/2y/w), where n = 0,1, 2,.... Some of the Hermite polynomials of 
low-order are listed in Table I. 

Now we are left with only the last terms of Equation (IV.33) which will give 
a solution for P(z). According to Equation (IV.40), we know the fact that 


4m 


Ui 




u 2 


4 n 
w 2 


By using these two relations, we rewrite the r-independent terms as 

9 jh 8P 4 

- 2k-—~ = Ul + u 2 = —(m + n). (IV.42) 

q oz w 2 
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H 0 {s) = 1 

i?i(s) = 2s 

H 2 (s ) = As 2 -2 

H 3 (s ) = 8s 3 - 12s 

H 4 (s) = 16s 4 - 48s 2 + 12 


Table 1. Hermite polynomials. From [16] 


Substituting the middle terms of Equation (IV. 18) into l/q, we obtain 


dP i 2 , iz , zo 

TP = - - ip m + n) = -p-y - (m + n + 1 -p;- 

q kw z z 2 + Zq z z z + z 0 2 

where 2/{kw 2 ) = z 0 /(z 2 + y 0 2 )- After integrating, we obtain 


iP = In l - 1 * 


i(m + n + 1) tan 


(IV.43) 


(IV.44) 


w / z 0 

which is similar to Equation (IV.22) for the fundamental mode. Combining all 
sources, such as Equation (IV.41) for g, Equation (IV.44) for iP and Equation (IV.35) 
for l/q, we finally have the following solution for higher-order modes: 


/ - r HlTT 

Ym,n rn 

W 


2-1 H„ 

Ws 


2-) e~ r2/v 
w 


exp 


Ay/ ✓ v _1/6 

—— — [m + n + 1) tan — 
2 R zq 


Utilizing this solution, Equation (IV.4) can be written in the form 

, < 

2 —i e -” 3 /'" 2 


_ WOrr 
^ m.n flm 

W 


2— ) Hr, 

Ws 


w 


exp 


Amplitude 


(IV.45) 


\ 


kr 2 


kz + —zr — (m + n + 1) tan — 
2it Zq 

L V Phase (ft ) 


(IV.46) 

where the positive integers m and n identify the modes. These integer indices deter¬ 
mine the shape of the beam intensity in the x and y direction, respectively. Several 
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examples of higher-order modes are shown in Figure 15. The intensity maxima are 
separated by anti-nodal lines. Note that mode label m equals the number of intensity 
anti-nodes proceeding from left to right, and n the number proceeding from top to 
bottom. 



-►X « 01 Wj 2 W 32 

Figure 15. Intensity patterns for Hermite Gaussian modes in the transverse plane, x 
and y. 
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V. FEL THEORY WITH HERMITE GAUSSIAN 

MODES 


In the previous chapter, we found solutions of a wave equation for a Hermite 
Gaussian (HG) optical beam in free space. In the actual FEL operation, a beam of 
electrons interacts with the Gaussian optical beam, which is bounded by the spherical 
resonator mirrors, as shown in Figure 16. 



Figure 16. A schematic of the FEL oscillator shows the Gaussian beam (blue), with 
radius w(z) expanding away from a waist wo centered in the undulator (green) of 
length L. Also, the optical phase o(z) (black) changes along L. One electron path, 
on-axis at r = 0, is shown, and all transverse dimensions are exaggerated. After [24] 


This chapter explores the FEL interaction, considering the characteristics of HG res¬ 
onator modes. We will find a modified pendulum equation describing the microscopic 
electron evolution and determine the design of the HG mode which can maximize the 
interaction with the electron beam. 
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A. MODIFIED PENDULUM EQUATION 

The electric field at the mode center (r = 0) in a fundamental Hermite Gaus¬ 
sian mode as derived in Chapter IV is 


E{ 0, z) = E 0 ^~ exp [i (kz + 4>{z ))}, 
w[z) 


(V.l) 


where E 0 is the electric field amplitude at the mode waist, k = 2tt/A is the mode wave 
number, and <p(z) = — tan -1 [{z — L/2)/z 0 \ is the mode’s phase. Also recall that the 
mode radius is w(z) = Wq[ 1 + (z — L /2) 2 / Zq 2 ^ 2 , where z 0 = 7rw 0 2 /A is the Rayleigh 
length, and wq is the mode waist radius at z = L/2. Although we will only deal 
with the fundamental Gaussian mode, results can easily be extended to higher-order 
modes [24], The mode is centered at r w = 1/2 (z = L/2), since this point is near a 
broad maximum in gain [24], Using dimensionless time t — z/L — ct/L in w(z) and 
4>{z) gives 


w 2 (t) = Wq 2 



(V.2) 


and 


<j)(r) = — tan 1 




(V.3) 


Recall the pendulum equation derived in the previous chapter: 


00 

C = v = “M COS [c + <A(r)], 


(V.4) 


where a(r) is the dimensionless field amplitude in Equation (V.l) by a(r ) = ci 0 w 0 /w(t), 
where a 0 = inNeKLEo/^rnc 2 is the dimensionless optical wave amplitude. Insert¬ 
ing Equation (V.3) into (V.4) gives us the modified pendulum equation, describing 
an electron’s evolution in an HG mode. The parameter L/z 0 compares the undulator 
length L with the Rayleigh length z 0 . The Gaussian beam radius w and phase 0 
depend on L/z$, as we saw in Equation (V.2) and (V.3). In order to illustrate some 
of the effects of the Gaussian beam, consider L/z 0 <C 1 so that we almost have plane 



waves [24], To lowest order in L/zo, we can then attain the following relations with 
L/zq < 1: 

UJ (t ) W 0 , (V.5) 


thus 


a(r) = ao 


_Wo_ 

w(t) 


ao, 


(V.6) 


and 



(V.7) 


using the small-angle approximation. Substituting Equations (V.6) and (V.7) into 
(V.4) yields the modified pendulum equation for the case of low gain (j < 7r) 


=u 


a 0 cos 




ao cos 
field 


(o + 1 k +T 


new phase new resonance 



(V.8) 


where C = Co + iy o T , the solution of C to zeroth-order in ao, has been used in weak 
optical fields (a 0 < 7r). The electrons have a new initial phase. But the shift Co — > Co + 
L/2z () in each phase is not significant since the beam is uniformly spread over each op¬ 
tical wavelength [24]. The gain spectrum has the form [2 — 2 cos u 0 * — z/ 0 *(sin n 0 *)] /n 0 * 3 
derived in Chapter III, where now = u 0 — L/z Q in an HG mode. The gain spectrum 

is now shifted away from the exact resonance by an amount of L/z 0 . 


B. GAIN IN GAUSSIAN MODES 

In a weak optical field ao < n and for low gain j < ir , the gain in the HG 
modes becomes 

G(t) oc jF g , (V.9) 

where j is the dimensionless current density, and Fq is the filling factor for HG modes. 
Now, let us find the maximum filling factor averaged over the undulator length L in 
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Consequently, the maximum filling factor, Fc max , occurs at z 0 max = L/y/l2 ~ (0.3)L. 

We have analytically derived that the gain spectrum is shifted by u 0 —> u 0 — 
L/ Zq when Lj z 0 -C 1. We have also shown analytically that the optimum gain occurs 
when L/ Zq ~ 3. Numerical simulations have shown that the shift in the gain spectrum 
still occurs for values of L/zq up to ~ 10, and is still roughly given by the approximate 
shift uq —► Uq — L/zq [26]. 
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VI. HYBRID (HG-WAVEGUIDE) RESONATOR 

MODES 


For short wavelength operation of an FEL, an open resonator bounded by 
spherical reflecting mirrors is feasible, since the laser beam maintains a small cross- 
sectional area as it evolves along the interior of the undulator. However, when the 
FEL operates at long wavelengths, such as in the sub-millimeter or THz region, the 
expansion of the laser beam due to diffraction can be a problem, as shown in Figure 
17. 







(d) (e) (f) 

Figure 17. Surface plots of two Gaussian beams propagating along the 2 -axis without 
waveguides. Compared to the evolution of a short wavelength beam (Figures (a), (b) 
and (c) - 3/irn), the evolution of a long wavelength beam (Figures (d), (e) and (f) - 
30/im) obviously has a larger amount of expansion. 


To accommodate lager diffraction at longer wavelengths, we need to have a larger 
undulator gap g, but the larger magnet spacing reduces the strength of the magnetic 
field according to B = 1.95B mag exp(— irg/Xo), where Aq is the undulator period [2]. 
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Consequently, the gain G will be decreased proportionally: G oc D~. A remedy is to 
use a waveguide structure, shown in Figure 18, to confine that the radiation. This 
chapter analyzes the Hermite Gaussian-Waveguide resonator modes, which will be 
designated as Hybrid resonator modes. The modes include characteristics of both 
Gaussian and waveguide modes. In order to understand the attributes of the waveg¬ 
uide, we solve the wave equation and explore the phase and group velocities of the 
laser pulse propagating through the waveguide. 



Figure 18. Geometry of a waveguide resonator. The origin is located at the center of 
the waveguide, with y = 0 at the bottom of the waveguide. Curved surfaces represent 
the mirrors. From [22] 


A. SOLUTION OF WAVE EQUATION 

The Hybrid resonator modes have been studied in detail by Elias [17], Degnan 
[18], and others. This section follows Elias’s analysis, which is the basis of the UCSB 
FEL. To find solutions of a wave equation, assume that the electrical conductivity of 
the walls in Figure 18 are infinite. First find solutions to the wave equation inside a 
perfect, infinitely-long, rectangular waveguide in terms of Gaussian waves propagating 
along the z-axis. These manipulations are similar to the derivations of Chapter IV; 
the difference is that we need to consider waveguide properties in the ^-direction. 


42 



















1. Fundamental Mode 

Each cartesian component of the electric (or magnetic) field must be satisfied 
by the Helmholz equation presented in Chapter IV: 


(V 2 + k 2 ) u(x, y, z) = 0, (VI. 1) 

where k = lu/c = 2tt/A, and a harmonic time variation, e~ lU}t , has been assumed. 
Note that u(x, y, z) can represent either a component of the electromagnetic field or 
of the vector potential of the light. We are interested in finding solutions describing 
Gaussian waves propagating along the z-axis without phase variations along the y-axis 
[17]. We express u(x,y,z), an electromagnetic wave traveling along the ^-direction 
within the waveguide, as 

u(x, y, z) = ip(x, z)e lknZ sin (nny/b ), (VI.2) 


where ip(x, z) is a slowly-varying complex function, e lknZ is a fast-varying carrier 
wave, sin(n7ry/&), b is the distance between the top and bottom plates, and n = 
±1, ±2, ±3,... is an integer in order to match boundary conditions at the walls of 
the waveguide (note y = 0 at the bottom wall). Inserting Equation (VI.2) into the 
scalar wave equation and dividing by e lknZ sm(mry/b) yields 


d 2 ip d 2 xp 


dx 2 


2ik ^ 

oz z oz 






= 0 . 


From this equation, we find that 


k n = \k 2 - 


n 2 n 2 


(VI.3) 


(VI-4) 


When we apply the paraxial approximation [12], d 2 'ip/dz 2 <C d 2 i\)/dx 2 , and we arrive 
at the paraxial wave equation 


d 2 y> 


~\~ < 2ik 1 


dip 


dx 2 ' n dz 


Let us assume a trial solution of the form 


= 0 . 


ip(x, z) = exp 


fW 


2 q(z) 


X 


(VI.5) 


(VI.6) 
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similar to the Hermite-Gaussian case, except that here there is no y variable. We need 
to find the unknown functions P(z) and q(z) in a similar method to that described 
in Chapter IV. After inserting Equation (VI.6) into (VI.5) and comparing terms of 
equal power in x, we obtain the following relations 


dq 

dz 


= 1 , 


(VI.7) 


(VI-8) 


and 

dP _ i 
dz 2 q 

In order to find P(z) and q(z), we use the definitions of w(z) and R(z) again, but 
these are not quite the same as the previous parameters since a new Rayleigh length, 
zr = w 0 2 k n /2 [17], differs from the previous one in Chapter IV, z 0 = 7rw 0 2 /A = 
wo 2 k/2. Using this new definition of zr, we integrate Equation (VI.7), yielding 


q = z + q 0 = z- i 


. n'o k 71 


= Z- IZr, 


(VI.9) 


where the complex beam parameter at the waist is qo = —iwo 2 k n /2. Starting with 
Equation (VI.9), q = z—izR , and employing Equations (IV.15) and (IV.16) in Chapter 
IV, we obtain the following relation 


2 q 


2 R + w 2 ' 


(VI.10) 


To determine the P(z), we insert Equation (VI.9) into (VI.8) and integrate the result, 
and we then have 


iP(z) = In 


w 


Wo 


n 

2 


-tan 


-l 


zr 


(VI.11) 


using the definition of w(z). Substituting Equations (VI.10) and (VI.11) into the trial 
solution, Equation (VI.6), we can finally get the solution of the fundamental mode 
inside the waveguide 


^{x,z) 





(VI.12) 
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This solution for the fundamental mode is also similar to the Equation (IV.24) in 
Chapter IV, but now k n appears instead of k, and the Gaussian form is only in the 
x-direction. Using this solution, Equation (VI.2) can be written in the form 




( \ 


/ \ W 0 —^. 2/^2 

u(x,y,z) = \ —e ' exp 

V W 

i 

, k n x 2 1 2 

k n z -\ -—-tan — 

2 R 2 z R 


Amplitude 


V 

\ Phase (t> ) 



tut 

sin(—?/)• (VI. 13) 


2. Higher-Order Modes 

The higher-order solutions in the waveguide also take the form of Hermite- 
Gaussian functions in Cartesian coordinates, as we now show. Assume a trial solution 
to Equation (VI.5) in the form 


ip — g exp 


P+^x 2 

2 q 


(VI. 14) 


where g is a function of x/w. We can easily find the derivatives of the unknown 
function g when we follow the notation mentioned in Chapter IV. They are 


dx 2 “ 


1 // . 2 / ik n \ t fik n \ fk n 

—g + - —x ) g + [ — ) g - — x g 

~ 2 w V q J V Q J \Q 


w 


and 


d ± = ^p+^p) 

dz 


1 ,dw .dP ik n x 2 dq' 
w > 9 lh X + 9% ~dz ~ 9 ^hf~dz) ’ 


(VI.15) 


(VI. 16) 


where (..)' = d(..)/d(x/w). Inserting Equations (VI.15) and (VI.16) into (VI.5) and 
dividing by g exp (iP + ik n x 2 /2q), we arrive at the following result 


1 g" 2 ik n xg' dw/dzg' 2 

0 I Ztrb n X ~ I X 

w z g wq g w z g 

^ _ ✓ V 


kn 


k 2 dq 


q 2 dz 


ik n 


dP 


— —7T + “-VT VT’ ■+- i —— — 2k n ~— = 0. 


dz 


x only = iii 


x 2 only = 0 


no x = u\ 


(VI. 17) 

Utilizing the same notation as in Chapter IV allows us to have g(x/w), P(z), and 
q{z) such that 

g(-)=H m (-), m = 0,1,2,..., (VI. 18) 

\wJ \wJ 
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iP = In 


in 

w ' 2 
w 0 


— i I m -\— I tan 1 —, 
2 / z 0 


and 


1 1 


2 i 


(VI.19) 


(VI.20) 


q R k n w 2 

Substituting these three equations into Equation (VI. 14), we finally have the following 
solution for higher-order modes 


lw o 


^Pm,n \ Um ( ) ® 

W \W 


X 


—x 2 /w 2 


exp 


k n x 2 ( 1\ z 

!| w( m+ V* 8 "' ^ 


(vi.21) 


Utilizing this solution, we can express the optical beam, Equation (VI.2), as 

/ \ 



H m (- 

w \w 


e -x 2 lw 2 exp 


Amplitude 


h -r 2 

k n z - 

n 2 R 


m + - 1 tan 


-1 


zr 


\ 


Phase 


. / TITT 

sin y—y 
(VI.22) 
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Figure 19. Intensity patterns E* mn for Hybrid modes in the transverse plane, x and 

y- 


46 


























Note that n, the waveguide constant in the sine function, differs from the previous 
n of H n in Chapter IV. However, this n together with m determines the shape of 
the mode, as shown in Figure 19. We chose the sm(nny/b) as a waveguide term in 
Equation (VI.2) to satisfy the boundary conditions; there are no fields for n = 0. 

3. Solutions for E n mn Modes 

Elias and Gallardo [17] proposed two sets of linearly independent solutions 
differing by the choice of sm(niry/b ) or cos(ran//6). In the solutions for sm(nny/b), 
which will be designated as E n mn modes, the dominant electric fields are parallel to 
the wide dimension of the waveguide, as shown in Figure 20. 



Figure 20. Electric field lines at the waist of the waveguide, which shows the E w 0 i 
mode. 

The solutions for the choice cos (mn//6), which will be designated as E a mn modes, 
the dominant electric fields are those which are perpendicular to the ^-direction of 
the waveguide not shown in Figure 20. When we chose the .E' 7r mn modes, waveguide 
losses are negligible since the electric fields near the plates are close to zero. In an 
FEL, the undulator magnets would be placed along the y-axis, so that an electron 
wiggles along the x-axis and couples to the electric fields of the mode, as shown in 
Figure 21. 
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Figure 21. Schematic of a Hybrid resonator, showing the location of the critical com¬ 
ponents after [19]. The hybrid waveguide consists of a set of two parallel conducting 
plates separated by a distance b. A relativistic electron wiggles along the x-axis. A 
hole may be placed in the middle of the right mirror to allow output of the laser beam 
(KAERI is using this type of mirror for their FIR FEL [20]). Otherwise, a partially 
transparent mirror will be used for the out-coupling. In case of entry and exit of 
the electron beam, a chicane is usually used to bend the electron beam around the 
mirrors. 


We are now ready to explore the solutions for E n mn modes. Recall the pro¬ 
posed solution u(x, y, z ) = ^(x, z ) sm(mry/b)e iknZ . We can then express the dominant 
electric fields E x as 

E x « AoVw sin (VI.23) 

where A 0 is a constant. Using this equation and Maxwell’s equation V • E = 0 in free 
space, we can obtain the E z component: 


E z 


A a d'ip mn 
ik n dx 


sin 



g ik„z 


(VI.24) 


Note that the E x component is much greater than E y and E z , i.e., E x E y , E z . Also, 
E y = 0, as represented by Figure 20. In order to determine the magnetic fields, we 
use Maxwell’s equation V x E = —dH/cdt. When we consider only the y component 
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of this equation, we have 


dBy 

cdt 


0E X 

dz 


-(ik n )A 0 i> mn sin 



ik n Z 
c ? 


and 

By « ^-A^ mn sin (^y) e ik " z , (VI.25) 

where dd> rnn /dz has been ignored, since the amplitude of the longitudinal components 
of the fields are much smaller than the corresponding transverse components, and the 
harmonic variation e~ %ut has been assumed, respectively. Next, let us discuss the 2 
component of V x E = — <9B /cdt. We then have 


dB z _ dE x 
cdt dy 


mi 

~b 


A 0 r/w cos 



ik n z 


and 

B z ~ cos e ik " z , (VI.26) 

where a harmonic variation e~ l0Jt has been assumed again. Lastly, the boundary 
conditions applied on the E^ mn modes are E x , E z = 0 at y = 0 and y = b. A 
summary of the E n mn modes solutions is given by the following table. 


E x « A 0 xp mn sm(nny/b)e iknZ 
Ey = 0 

E z « - (. A 0 /ik n ) (d'lpmn/dx) sm(nny/b)e tknZ 


Table 2. E n mn - modes solutions 


Starting with these solutions, we find a relation between the coefficient A 0 and 
the power of the electromagnetic field. The expression for total energy density [21] 
is: 

u=^(E 2 + B 2 ). (VI.27) 
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B. PHASE AND GROUP VELOCITY 

The general expressions for phase and group velocities of a electromagnetic 
wave are v p = u/k, and v g = d u/ dfc.Then from Equation (VI.4) with z <C zr, we 
have: 

k n 2 1 1 /nn\ 2 

Ll > 2 Vpn 2 C 2 Vcjfe/ 

Therefore, 

v pn = - C (VI.32) 

y 1 — ( nnc/cub ) 2 

Note that the phase velocity of the wave depends on the waveguide number, n = 
1,2,3,..., the frequency, ui, and the height of the waveguide, b. Also, the phase 
velocity, v pn , is faster than the speed of light, c, since 
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The group velocity of the wave can be derived from Equation (VI.4). Differentiating 
this equation with respect to u, we get 


2 k r 


d k n 
den 


UJ 


= 2 — 


dk n 


LU 


den k n c 2 


Hence, the group velocity is 


den klc c 2 


Vgn dkn 


en v, 


pn 


We can rewrite Equations (VI.43) and (VI.45) as 
c c 


Vpn 


1 — ( nnc/uby 


1 — (cnX/ZXfby 


1 - (n\/2by 


and 


Vgn ~~ ~ ~ C ^ 1 l 2b) ' 


J pn 


Thus, v pn v gn = (? is also satisfied inside the waveguide. 


(VI.33) 


(VI.34) 


(VI.35) 


(VI.36) 


Note that when the wavelength is equal to the cut-off wavelength, i.e., A c = 2b 
and 9 = 90° as shown in Figure 22, the phase velocity is infinite and the group velocity 
is zero so that no energy propagates along the waveguide. 



Figure 22. A sinusoidal electromagnetic wave entering one end of the waveguide with 
its electric field parallel to the wide dimension and propagating at an angle 9 will be 
successively reflected by the top and bottom walls and follow a zig-zag path. A is 
near the cut-off wavelength. From [23] 
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The cut-off frequency f c is given by 
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VII. 


FEL THEORY WITH HYBRID 
(HG- WAVEGUIDE) MODES 


A Hybrid mode occurs in a waveguide consisting of two flat conducting plates 
seperated by a distance b, each parallel to the x-z plane. This waveguide serves as 
an optical beam pipe between the resonator mirrors. In principle, there are image 
charges in the conducting walls that may affect the electron trajectories, but they 
are very small. More important things are the space charge forces in the beam itself 
which are assumed negligible here, as in most FELs. In this chapter, we will explore 
a modified pendulum equation appropriate to Hybrid modes. We will then design the 
relevant FEL which results in maximum energy gain, using the same analysis as in 
Chapter V. 


A. MODIFIED PENDULUM EQUATION 


The electric field at the mode center (x = 0, y = 6/2) in a fundamental Hybrid 
mode as derived in Chapter VI is 

E (o, | ¥ z) = EoJ^exp [i {k x z + ^>{z ))\, (VII. 1) 

where E 0 is the electric field amplitude at the mode waist, k\ = ^Jk 2 — n 2 /b 2 , and 
(f){z) = — 0.5tan -1 \{z — L/2)/zr], where zr = wo 2 ki/2 is the Rayleigh length for 
this Hybrid mode case, and w 0 is the mode waist radius at z = L/2. Although we 
will only deal with the fundamental Hybrid mode, results can easily be extended to 
higher-order modes. We can again introduce the dimensionless time, r = z/L, so 
that the mode is centered at t w = 1/2, near a broad maximum in gain [24], Through 
a derivation similar to that in Chapter V, we have 

i+ (S’H)} (viL2) 

and 



0(r) 




(VII.3) 
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Our original pendulum equation in Chapter III is 

00 ~ 

C = v = o(t) cos [( + <j>(r)\ , (VII.4) 

where a(r) is the amplitude factor in Equation (VII. 1), where a 0 = AirNeK LEq/ ^mc 2 
is the dimensionless optical wave amplitude. Using the approximation, L/zj; C 1 de¬ 
scribing a near plane wave, gives us to lowest order in L/zr 

w{t ) ~ wo, (VII. 5 ) 


so that 


and 


a(r) = cio 


wo 
w(t ) 


a 0 , 


0(r) « - 


L 


2z 


R 


1 

T ~2 


(VII.6) 


(VII. 7) 


Substituting Equation (VII.6) and (VII.7) into (VII.4) yields the modified pendulum 
equation for the case of low gain (j < 7r), derived here for the first time: 


^ =u m a 0 cos 


Co + "or - 


4 Zr 


T — 


a 0 cos 
field 


Co 


4z 


L ( L 

+ T I/ 0 - 


R 


2 Zr 


new phase new resonance 


(VII.8) 


where C = Co + l/ i)T has been used assuming weak optical fields (ao < 7 r). The electron 
has a new initial phase C,o + L /2zr, but this shift can be ignored because the electrons 
are uniformly spread along each optical wavelength. However, we cannot ignore the 
new resonance parameter, z/ 0 —> (uq — L/2zr), which causes a measurable change in 
the FEL operation. In a weak optical field a 0 < it, the gain spectrum has the form 
[2 — 2 cos vq* ~ ^o*(sin z/ 0 *)] /n 0 * 3 derived in Chapter III, where now v 0 * = v 0 — L/2 zr 
in the Hybrid mode. 
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B. GAIN IN HYBRID MODES 


In a fundamental Hybrid mode, the transverse mode area can be used to 
calculate the filling factor. This affects the gain, since gain is proportional to the 
filling factor (G oc jF H ). The cross-section area of the elliptical Hybrid mode can be 
expressed nr x r y , where r x = w(z), r y ~ 6 /2, and 6 is the height of the waveguide as 
shown in Figure 23. 




Figure 23. Cross-section of the fundamental Hybrid mode. The area of the elliptical 
cross section of the mode is irr x r y ~ nw(z)b/2. 


By using a similar analysis as in Chapter V, we can find what zr value gives 
the maximum Fr, which leads to the maximum gain. In a weak optical field a 0 < 7r, 
the gain for the fundamental Hybrid mode is given by 

G(r) oc jF„, (VII.9) 


where j is the dimensionless current density, and Fr is the new filling factor for the 
Hybrid mode, which we will explore in detail through this section. The averaged 
filling factor Fr becomes 


electron beam area 


27 m 2 


Fr = —— — —— = (VII. 10) 

optical beam Area irbw(z) 

where is the electron beam radius (assumed constant), w(z) is the optical mode 

semi-major axis along x, and 6/2 is the optical mode semi-minor axis along y. We 

can write the area of the fundamental Hybrid mode as 


... 7t6 , . 7t6 / 

A(z) = —w(z) = —w 0 1 -|-- 

Zr 2 


1/2 


(VII.ll) 
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where zr = w 0 2 ki/2, and k\ = ( k 2 — n 2 /b 2 ) 1 ^ 2 . Thus, the averaged mode area be¬ 
comes 



nbwQZR 

2 L 


7 rb 

VWl 


A{z) dz = 



2 z R \dz R 2 

Vz^Lf L 2 




(VII. 12) 


where vjq = (2zrJ ki ) } ' 2 has been used in the last step. We need to find the ZR max 
value which leads to the maximum filling factor, and minimum A. However, d A/ d zr 
is a complicated, so let us find ZR max numerically. From the numerical solution, we 
obtain the estimated ZR max value, which yields the smallest optical area, as shown in 
Figure 24. 



Figure 24. Averaged optical mode area A for fundamental Hybrid mode versus 
Rayleigh length zr, when b = 1cm and L = lm. We see from this graph that 
F H max occurs at z# max ~ (0.4)L. 


56 




















Therefore, the maximum filling factor, F H max , occurs at zu " ax « (0.4)L, giving a peak 
value of G/jF H (gain normalized to j and Fh) of about 0.135 and at z/ 0 rnax ~ 3.85. 
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VIII. 


CONCLUSION 


In this thesis, we have reviewed the basic FEL theory and found solutions of 
the wave equation for Hermite Gaussian (HG) resonator modes in free space. We 
have then analyzed the effects of HG modes on FEL theory, including the pendulum 
equation and gain. We have then studied the Hybrid (HG-Waveguide) resonator 
modes. There are crucial for long FEL wavelength, where a waveguide must be used 
in the resonator cavity to prevent the optical beam from diffracting into contact with 
the undulator. We have analyzed, for the first time, FEL operation with these Hybrid 
modes. 

A waveguide in the resonator cavity produces a change in the FEL resonance 
parameter. This change has the effect of shifting the gain spectrum away from the 
traditional (free space) value of resonance. In a fundamental Hybrid mode, the trans¬ 
verse mode area must be used to compute the filling factor, which also affects the 
gain. We have found a new optimized Rayleigh length, associated with the waveguide, 
which allows the FEL to have maximum gain. Finally, we recommend the following 
future work: 

• Examine the beam coupling to free space mode outside the Hybrid mode 
region. 

• Examine the output coupling with a partially transmissive mirror versus beam 
output that using a mirror with a hole. 
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